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Abstract. A linear but global numerical model for the magnetorotational instability (MRI) in disk geometry is 
considered to estimate the instability parameters for galaxies. Similarity rules suggested by a local analysis are 
applied to reveal a universal behavior approached by the results of global calculations for large magnetic Prandtl 
numbers. The findings are used to estimate the MRI characteristics for galaxies with their very large magnetic 
Prandtl numbers which cannot be attained in any numerical simulations. The resulting minimum field for the 



instability, B n 



10" 



G, is small compared to any seed fields currently discussed. The growth times of MRI 



are estimated to be on the order of the rotation period of the inner rigidly-rotating core, i.e. 
MRI excites preferentially the magnetic field modes of quadrupolar symmetry. 
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1. Introduction 

The magnetorotational instability (MRI) has been discov- 
ered by Velikhov (1959) for hydromagnetic Couette flows; 
in the astrophysical context it has been established by 
Fricke (1969) and Balbus & Hawley (1991). Several prop- 
erties of the instability make it also relevant for the theory 
of galactic magnetism. 

The instability exists due to magnetic fields, but the 
minimum field for MRI is really very small. We shall see 
that for galactic conditions the minimum field, -B m in — 
10 -25 G, is much smaller than the seed magnetic fields of 
protogalaxies whatever origin of the seed fields we assume. 
With its growth time of order the rotation time the MRI is 
fast, with e-folding times below 0.1 Gyr. The observation 
of Faraday rotation of high-redshift objects by Athreya et 
al. (1998) and Carilli & Taylor (2002) seems to indicate the 
real existence of such short growth times of fields (Lesch 
& Hanasz 2003). 

It was recognized by Beck et al. (1994) that the seed 
field can be largely amplified by interstellar turbulence 
before the global hydromagnetic dynamo starts operating 
(Cattaneo 1999, Schekochihin et al. 2004). Turbulent ve- 
locities of about 10 km/s are indeed indicated by spectral 
line broadening of protogalactic disks (Wolfe 1988). The 
origin of the turbulence is, however, uncertain. Supernova 
explosions would drive a much weaker turbulence (Ferriere 
1992, Kaisig et al. 1993). Another driver is MRI which was 



demonstrated by Sellwood & Balbus (1999) to produce its 
own turbulence in the nonlinear regime. 

A more specific problem of galactic dynamos to which 
MRI can also be relevant is the parity selection by the 
global magnetic dynamo. Galactic dynamo models prefer 
the global fields symmetric about the disk midplane (Stix 
1975, and many other references in Riidiger & Hollerbach 
2004) which symmetry type is favored by observations 
also. In order to be amplified, however, the quadrupolar 
field should already be present in the seed field for the 
dynamo. Not all origins for the seed field currently dis- 
cussed lead to the required symmetry, e.g., battery effect 
of Biermann supplies dipolar global fields (Krause & Beck 
1998). We shall demonstrate that MRI excites preferen- 
tially global modes with quadrupolar parity even if the 
(axial) destabilizing field is antisymmetric with respect to 
the disk midplane which property suggests MRI as a nat- 
ural supplier of the seed field for global galactic dynamos. 

One could argue that MRI in the weak-field regime 
is rather local by nature and not efficient in producing 
global patterns. The local theory of ideal fluids yields 
the growth rate, 7, for MRI proportional to the Alfven 
frequency, 7 ~ VaA in the weak-field case, Va <C fl/k 
(Va = B /y/[i p is the Alfven velocity, k is the vertical 
wave number, Q is the local angular velocity, and Bq is 
the background axial magnetic field). The growth rate 
seemingly increases with the wave number. The growth 
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is, however, stabilized by the microscopic magnetic diffu- 
sivity, 77, so that 7 4- r\k 2 ~ V^k. Marginal stability then 
demands Va — T]k. The smallest field strength produc- 
ing the instability belongs to the smallest possible k. The 
low bound on magnetic field producing MRI is therefore 
controlled by the global geometry. 

We try to find out in the present paper whether MRI 
should be expected for protogalaxies where only very weak 
seed fields are present and the magnetic Prandtl number 
is extremely large (see Kulsrud & Anderson 1992). The 
extent to which the seed field can be amplified by the 
turbulence is evaluated. The preferred geometry of mag- 
netic field in the global unstable modes is determined. To 
this end the disk model of Kitchatinov & Mazur (1997) 
is used which is outlined in the next Section. Section 
presents the results of the computations. The implications 
for galactic dynamos are discussed in the final Section 0] 

2. The model 

The model concerns a rotating disk of constant thickness, 
2H, threaded by a uniform axial magnetic field. The rota- 
tion axis is normal to the disk and the angular velocity, Q, 
depends on the distance, s, to the axis. This dependence 
is parameterized by 



with 
O(s) = 



-l/r. 



(1) 



(2) 



This profile describes almost uniform rotation up to the 
distance sq, and the law £1 ~ £Io(sq/s), implying constant 
rotation velocity for distances s sq, as is observed for 
nearly all galaxies. We use n = 2 in Eq. and so/H = 5 
for the aspect ratio. 

We assume the fluid be incompressible, div U = 0. 
The fluid compressibility is probably not important for 
our analysis because MRI proceeds via noncompressive 
Alfven-type disturbances. The magnetic field and the flow 
in the disk obey the induction equation, 



dB 

~dt 



curl (U x B) + i]AB, 



and the equation for vorticity, W = curl U, 

at 



= curl (U x W + J x B/p) + i/AW, 



(3) 



(4) 



which results from taking the curl of the momentum equa- 
tion to eliminate the pressure. J — (l///o) curl_B is the 
current density. 

To perform the linear stability analysis, we linearize 
J2J and (0J about the rotation and the uniform axial 
field, Bo = Bq€. z (e z is the unit vector along the rotation 
axis), and introduce the normalized variables, 

b = B/B , j = ^-J, u = U/(Hn ), u = w/n (5) 
-on 



This leads to the three basic dimensionless parameters of 
magnetic Reynolds number, Rm, Hartmann number, Ha, 
and magnetic Prandtl number, Pm, 

n H 2 b q h v a h v 

Rm = , Ha = — ^=^= = , Pm = — . (6) 

They control the equation system for the normalized dis- 
turbances, 



Rm curl 



(s£l(s) 



x b — e z x u ] + Ab. 

2 



db 

~di 

= Rm curl ( sCl(s) x u) = e z x u 

dt \ v ; 20(s) 



Pm Ha 2 
Rm 



curl (e z x j) + PmAw, 



where k is the normalized epicycle frequency 
, 2Cl ds 2 £l 

K = A—' 

s as 



(7) 



(8) 



time and distances are normalized to the diffusive time, 
H 2 /rj, and the half-thickness, H, of the disk. The bound- 
ary conditions on the disk surfaces are (i) stress-free for 
the flow and (ii) pseudo-vacuum condition for the mag- 
netic disturbances, i.e. 



x b = 0. 



(9) 



The solutions are required to be regular at the rotation 
axis and to vanish at infinity. 

The linear stability analysis with b, oj ~ cxp(7i) leads 
to an eigenvalue problem for the equation system (JJJ) . The 
problem has been solved numerically with the method de- 
veloped by Kitchatinov & Mazur (1997). As usual, the 
system (J7J) allows two types of solutions with different 
symmetries about the disk midplane. One of the sym- 
metry types combines a symmetric magnetic field with 
an antisymmetric flow field. This symmetry type is called 
"quadrupolar" (we are mainly interested in the symme- 
try of the magnetic pattern). The other symmetry type 
of "dipolar" perturbations combines antisymmetric mag- 
netic fields with symmetric flows. The eigenmodes are also 
distinguished by their azimuthal wave number, to, because 
the axisymmetric modes with to = and nonaxisymmet- 
ric modes with to > 1 satisfy independent equation sys- 
tems. 

The primary aim of the linear theory is to define the 
stability boundary in parameter space which separates the 
region of stable perturbations with negative or zero real 
part of 7 from the instability region with (exponentially) 
growing perturbations. The stability boundary depends on 
the symmetry type of the excitations. The growth rates 
of the unstable disturbances are less significant because 
they characterize those states of the system which have 
little chance to exist just because they are unstable. The 
instability region is the realm of nonlinear simulations. 
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Fig. 1. The stability diagram for axisymmetric modes of 
quadrupole parity. The right branches and minima of the 
neutral stability lines remain constant for sufficiently large 
magnetic Prandtl numbers. 



3. Results 

Because of the galactic implications, we focus on the model 
results for large magnetic Prandtl numbers, Pm 3> 1. 
Figure ^ shows the neutral stability lines for axisymmet- 
ric quadrupolar modes. As discussed below, modes of this 
symmetry type are excited preferentially. The instability 
region is placed between the left and right branches of the 
bounding lines which branches define the minimum and 
maximum field strengths for the instability, respectively. 
The minima of the lines define the minimum Reynolds 
number for the instability. 

Several neutral stability lines for Pm increasing from 
1 to 1000 are shown in Fig.^ Our numerical model can- 
not match the huge values of Pm expected for galaxies. 
However, Eqs. (|A.6Jl - (|A.8|l suggest (and the results of 
Fig. H confirm) that the neutral stability lines approach 
a universal behavior in the high Pm limit. The behavior 
obeys the so-called similarity rules from which the MRI 
characteristics for arbitrarily large Pm can be estimated. 

Three similarity rules can be found from the numerical 
results illustrated by Fig. Q 

(i) On the right-hand branch the lines approach a con- 
stant slope for high Rm, so that one always has 



Rm 



n H 



2.4. 



(10) 



Obviously, with ~ c s one finds the magnetic Mach 

number slightly below unity. 



Fig. 2. Variation of the normalized growth rates 7T rot 
from the left-hand branch to the right-hand branch across 
the stability diagram of Fig. ^ for Rm/VPm = 100 {left) 
and Rm/^/Pm = 1000 (right). 



(ii) As Pm grows, the lines approach a common mini- 
mum of 



Rm 



(11) 



(iii) The left-hand branches on Fig. are close to ver- 
tical lines at large Rm. The lines are almost equidistant in 
the logarithmic scale suggesting that the Lundquist num- 
ber S is approaching a constant, i.e. 



S = 



V A H 



1.7. 



(12) 



The values l|10|l - H12|l provided by the global model do 
not coincide with the estimates (IA.6I) - I)A.8|) from a local 
treatment. The difference should be partly attributed to 
the different definition of the basic parameters in the local 
analysis (cf. eqs. © and HA.3l) b However, the scaling used 
for large Pm was inferred from the local analysis presented 
in the Appendix. 

Figure shows how the growth rates of unsta- 
ble perturbations vary along a horizontal line (constant 
Rm/VPni) from the left to the right branches of the sta- 
bility diagram. The growth rates are shown in units of the 
rotation frequency, t~ \ = £l /27r, of the disk core. It is 
suggested that as the Reynolds number increases the lines 
join together in a common dependence for different Pm. 
The typical growth times are of the order r rot . 

Figure |21 illustrates the structure of the global modes. 
The broken line encircles the rigidly rotating core, sq = 
5H. The modes are truly global with large horizontal 
scales. The scale, however, decreases with increasing am- 
plitude of Bq. This behavior is accompanied by an increase 
of the pitch angle between the azimuthal direction and 
the field vector. The pitch angle is small for the minimum 
magnetic field. It increases to almost 90° for the maximum 
field producing MRI. In both limits the magnetic-induced 
angular momentum transport vanishes. 
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Fig. 3. Magnetic field in the midplane of the disk for the quadrupole-parity mode at Rm/vPm = 1000. From the left 
panel to the right one: the marginally stable mode for the minimum magnetic field, the mode with maximum growth 
rate 7 = 1.03/T rot at Ha ~ 250, and the marginally stable mode for maximum magnetic field. The broken circle is the 
radius so = 5H of the rigidly rotating core of the disk. The background rotation is anticlockwise. Note the increase of 
the pitch angle and the decrease of the horizontal scale with increasing field strength. 



The linear theory cannot yield the magnitude of the 
magnetic stress but it can compare it with the kinetic 
stress. The angular momentum flux by Reynolds stress 
is also outward but it is small compared to the Maxwell 
stress (cf. Brandenburg et al. 1995; Riidiger et al. 1999). 
The relation reflects the magnetic nature of MRI. Also the 
energy of the unstable modes is dominated by its magnetic 
part. 

Figure 01 shows the stability diagram for the axisym- 
metric dipolar modes. Comparison with Fig. Q] reveals a 
preference of the quadrupole parity. The neutral stability 
lines of Fig.QJlay totally inside the correspondent instabil- 
ity regions of Fig. ^ This means that when the system ap- 
proaches instability to the disturbances of dipolar parity, 
it is already unstable to quadrupolar perturbations. In the 
same sense, the axisymmetric disturbances are preferred 
compared to the nonaxisymmetric ones (Kitchatinov & 
Mazur 1997). 



4. Implications for galactic dynamos 

4.1. Basic parameters 

To probe the relevance of the MRI for galactic magnetism, 
the parameters controlling the instability should be esti- 
mated for galaxies. 

We assume the angular velocity of the disk core as 
l^o — 3 • 10~ 15 s _1 which corresponds to a rotation period 
r rot ~ 70 Myr (Carignan et al. 1990; Sofue 1996). The 
half-thickness of the disk is assumed as H ~ 10 21 cm ~ 
300 pc. For a plasma with not too low degree of ionization, 
conductivity is controlled by electron scattering on ions. 
The magnetic diffusivity then is, 
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Fig. 4. The stability diagram for dipolar modes. The 
modes are harder to excite than the quadrupolar modes 
of Fig. 1. 



where In A ~ 20 is the Coulomb logarithm (Spitzer 1962). 
Accepting T ~ 10 4 K for the characteristic temperature, 
we find 



Tj ~ 10 7 cm 2 /s, 



Rm ~ 3 • 10^ 



(14) 



Viscosity is more difficult to estimate. The expression for 
the fully-ionized hydrogen (Spitzer 1962), 



10 13 T~ 3 / 2 



In A 
"20" 



cm 2 /s, 



(13) z/ = 6.5-10' 



T 5/2 



20 \ 2/ 

LTa ) cm /s ' 



(15) 
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can be applied to Hn-regions. With the proton density, 
n p ~ 1 cm -3 , it yields 

v ~ 7 • 10 17 cm 2 /s, Pm ~ 7 • 10 10 . (16) 

In the warm interstellar medium, however, neutrals can 
provide a larger viscosity, v n ~ vt^im ( w t is the thermal 
velocity and £ nn is the mean free path for neutral colli- 
sions). With the neutrals density, n n ~ 1 cm -3 , we find 

v ~ 10 21 cm 2 /s, Pm ~ 10 14 . (17) 

Fortunately, accuracy in estimating the viscosity is not im- 
portant. The viscosity value is needed to ensure only that 
the minimum ratio 1(11(1 for the instability is exceeded. On 
combining l|14(l with 1(16(1 or (|1 T|) , we confirm that the mag- 
netic Reynolds number exceeds the critical value resulting 
from (|ll|l by many orders of magnitude. This suggests 
that the range of magnetic field strengths producing MRI 
should be very broad. 

The maximum magnetic field for the instability results 
from ifTTijl . 

B max ~ 6 /xG, (18) 

(cf. Sellwood & Balbus 1999) which is surprisingly close 
to the observed strengths of galactic fields. We applied the 
estimation of Alfven velocity, Va — 2 B km/s with B in 
fiG, to obtain (fTHf) . 

The minimum field strength for MRI can be estimated 
from 1)12(1 to the extremely small value 

B min ~ 1(T 25 G. (19) 
4.2. Seed field amplification 

Magnetic field amplification by a hydromagnetic dynamo 
needs a seed field to start from. The origin of such a seed 
field for galaxies is very uncertain. Three possibilities can 
currently be found in the literature. First, the battery ef- 
fect of Biermann can produce the seed field. The corre- 
sponding amplitude was estimated by Davis & Widrow 
(2000) and Gnedin et al. (2000) to be about 10~ 20 G in col- 
lapsing protogalaxies. Second, collisional friction between 
ions and neutrals on ionization fronts may produce fields 
of order 10~ 16 G (Birk et al. 2002). Third, the cosmolog- 
ical magnetic field might be as strong as 10~ 10 G (Ratra 
1992). All the quoted field strengths can be further am- 
plified by a factor ~ 10 2 ...10 3 upon the collapsing phase. 

Whatever the origin of the seed field might be, the 
newly formed galaxies are well inside the instability re- 
gion between the limits given in Eqs. I|18|) and i|19[l . The 
region should be entered somehow, hence, the instabil- 
ity should start even before the galaxies are formed. MRI 
in young galaxies is operating in far supercritical region, 
B B m i n . Growth rates of order r rot are expected for this 
case. The maximum growth rates initially belong to small 
spatial scales, £<^H, with Q Q £ ~ Va » rj/£ (Fricke 1969; 
Balbus & Hawley 1991). The scale range certainly exists 
for the very large values given by 1)15(1. MRI at its non- 
linear stage is known to be able to amplify the magnetic 



field producing the instability. As the unstable perturba- 
tions grow, the increase of their magnetic field should shift 
the maximum growth rates to larger and larger scales until 
the global perturbations start dominating. 

The turbulent amplification of magnetic fields was 
studied by Beck et al. (1994). They did not specify the 
origin of turbulence in favor of quoting observational ev- 
idence for the turbulence. We argue here that the MHD 
turbulence can be produced by MRI. 

If we assume that the nonlinear amplification of mag- 
netic fields proceeds with about the same rate as the lin- 
ear growth, the seed field can be amplified by a factor of 
iamp(T-gai) ^ exp (r gal /r ro t ) for a galactic age r ga i. With 
the value flo ~ 3 ■ 10~ 15 s _1 it takes less than 2 Gyr 
to produce a 1 /iG field from even the weakest seed field 
expected from the Biermann battery effect. This may ex- 
plain regular magnetic field of fiG strength for galaxies 
with high redshift z ~ 2 (r ga ] ~ 1 ... 3 Gyr) . Restrictions 
on the magnitude of the seed field, however, would emerge 
if the microgauss fields were already present in very young 
galaxies as found at z ~ 3.4 by White et al. (1993). 

4.3. Field geometry 

We have seen in Section 01 that MRI produces prefer- 
entially the global perturbations of quadrupolar parity 
though the background field for the instability is dipolar. 
Therefore, MRI supplies the seed field of 'correct' parity 
for global dynamos. It may be mentioned also that the 
pitch angles (the angles between the azimuthal direction 
and the field vector) in aO-dynamo models are often too 
small. The small angles result from an efficient winding 
of magnetic field by nonuniform rotation. Figure [3] shows 
that the pitch angles for the most rapidly growing global 
modes of MRI are not small. 

It is characteristic for linear global MRI models that 
the symmetries of flow and field always differ. If the mag- 
netic field is symmetric with respect to the midplane 
then the flow is antisymmetric and v. v. Hence, if re- 
ally the global MRI is present in galaxies then the as- 
sociated flow field is antisymmetric with respect to the 
midplane. Then the effective rotation rate would differ 
between hemispheres. Nonlinear calculations must show 
whether in the final, saturated phase these characteristic 
mixed-mode properties of MRI-induced magnetism sur- 
vive. It might be important in this connection that in the 
evolving galaxy with SN-explosions a second driver of in- 
terstellar turbulence occurs which may finally dominate 
the magnetic-field production. 
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Appendix A: Local analysis of marginal stability 

The guidance for appropriate scaling of our numerical re- 
sults came from an elementary local analysis. The local 
treatment of MRI provides expectations and better un- 
derstanding for global simulations. 

We use a Cartesian coordinate system corotating with 
the local angular velocity, f2, and axes x,y and z point- 
ing in the radial, azimuthal and vertical directions re- 
spectively. The local approximation concerns perturba- 
tions with scales small compared to the global scale. Then, 
the rotation law can be approximated by the shear flow 
Uq = -e.yQ.qx where q is the (constant) local shear. The 
linearized MHD equations in the local approximation read 



qflB' x e y - rjAB' = 0, 



dB' du' n dB' 

-57 B — x 1 n ^- 

ot oz ay 

— — + 2\le z x u — xqll— q\lu, r e u 

dt H dy H x v 



B dB' 1 



VP' - vAu' = 0, 



(A.l) 



HoP dz p 

(Balbus & Hawley 1991; Brandenburg et al. 1995) where 
P' is the pressure fluctuation including the magnetic term, 
Bq is the amplitude of the background axial magnetic 
field, B = B e z . 

Considering plane waves with B',u',P' ~ exp(7i + 
ikz) leads to the dispersion equation, 

( 7 + V k 2 ) 2 ((7 + vk 2 f + 2 (2 - q) fi 2 ) + 

+ui% (ujI - 2qn 2 + 2 (7 + vk 2 ) (7 + rjk 2 )) = (A.2) 

with the Alfven frequency lua = kBo j JJhQp. 

We need the boundary separating the regions of sta- 
bility and instability in parameter space. The magnetic 
Reynolds number and Hartmann number have to be rede- 
fined for the local treatment as 

- l " " (A.3) 



Rm = 



rjk 2 



Ha = 



k' 2 . 



<r\v 



To find the marginal stability equation we put 7 = in 
Eq. (|A.2(1 . After some algebra it yields 



Rm 2 = 



Pm 2 fl + Ha 2 ) 2 



2 (q Pm Ha 2 - 2 + q) 



(A.4) 



With a fixed value of Pm, this equation leads to neutral 
stability lines similar to those of Fig. ^ For a given suf- 
ficiently large Rm, the instability exists in the region be- 
tween some minimum and maximum values of Ha. We call 
the lines defining the minimum and maximum Ha the left 
and the right branches of the marginal stability line. 
Then the minimum Reynolds number for instability is 



Rm 2 lin = Pm- [ 1 



at 



qPm 



qPm 



(A.5) 



(A.6) 



For large Pm, the ratio of Rm/VPm at minimum ap- 
proaches a constant and so does the correspondent Ha. 
This behavior at large Pm agrees with the numerical re- 
sults of Section|2| Actually, the scaling used in that Section 
was motivated by the local analysis. 

The right branch of the neutral stability lines at high 
enough Rm can be found by considering l|A.4|) for large 
Ha, 



Rm 2 = £^Ha 2 
2q 



(A.7) 



Again the ratio of Rm/vPm is an appropriate scaling for 
the slope of the right branch. 

The left branch is controlled by the diffusivities. Its 
expression for large Rm can be obtained by by setting the 
denominator in HA.4II to zero. 



Bi 



2;„2 



Hoprj k 



(A. 



The smallest magnetic field producing MRI belongs to the 
smallest possible wave number so that also the left branch 
is controlled by global disturbances. Equation (|A.8() shows 
that the minimum Lundquist number is constant on the 
left branch for sufficiently high Rm. 
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